Low-temperature thermodynamics of one-dimensional alternating-spin 

Heisenberg ferromagnets 
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Motivated by a novel bimetallic chain compound in which alternating magnetic centers are fer- 
romagnetically coupled, we investigate thermodynamic properties of one-dimensional spin-(S', s) 
Heisenberg ferromagnets both numerically and analytically. On the one hand, quantum Monte 
Carlo calculations illuminate the overall thermal behavior. The specific heat may exhibit a double- 
peaked structure at intermediate temperatures for S > 3s in general. On the other hand, a modified 
spin-wave theory precisely describes the low-temperature properties. Expanding the specific heat 
and the magnetic susceptibility, we reveal an analogy and a contrast between mixed-spin ferromag- 
nets and ferrimagnets. 

PACS numbers: 75.10.Jm, 75.40.Cx, 75.30.Ds, 75.40.Mg 



I. INTRODUCTION 

Much effort has been devoted to designing molecular 
systems ordering ferromagneticallyi One possible ap- 
proach consists of assembling molecular bricks so as to 
obtain a low-dimensional system with a nonzero resultant 
spin in the ground state and then coupling the chains or 
the layers again in a ferromagnetic fashion. Numerous 
heterospin chain compounds have been synthesized along 
this line. Gleizes and Verdaguer- made an attempt to al- 
ternate two types of metal ion along one crystallographic 
axis with antiferromagnetic intrachain interaction and 
obtained a pioneering example of quasi-one-dimensional 
ferrimagnets, of formula Mm4Cu(dto) 2 (H 2 0) 3 • 4.5H2O 
(A = Ni, Cu; dto = dithiooxalato = S 2 C 2 2 ). 
Kahn et al£ synthesized another series of bimetal- 
lic chain compounds ACu(pbaOH)(H 2 0) 3 
nH 2 (A - Fe, Co, Ni, Cu; pbaOH 
2-hydroxy-l, 3-propylenebis(oxamato) — CyH 6 N 2 07), 
one of which indeed attained the three-dimensional fer- 
romagnetic order at low temperatures^" Caneschi et al£ 
took an alternative strategy of bringing into interaction 
metal ions and stable organic radicals. This idea was de- 
veloped toward polymeric chain compounds^ The wide 
variety of chemical explorations stimulated the physical 
interest in mixed-spin chains! 7 ! 8 ! 9 ! 10 : 11 : 12 : 13 ! 1 ^ 15 ! 16 

Most of the thus-far synthesized heterospin systems 
are characterized as ferrimagnets. Ferromagnetic intra- 
chain coupling is observed in few cases. In such cir- 
cumstances, MnNi(N0 2 ) 4 (en) 2 (en = ethylenediamine = 
C 2 H 8 N 2 ) f i 7 . proved to be a quasi-one-dimcnsional mixed- 
spin ferromagnet 18 and caused renewed interest in mixed- 
spin chains. Gillon et alML calculated the spin density 
distribution by means of the density functional theory 
and quantitatively visualized the ferromagnetic nature 
of the Mn(II)-Ni(II) interaction. Fukushima et al~ per- 
formed high-temperature series expansion of the thermal 
quantities and argued the magnetic structure including 
single-ion anisotropy and interchain exchange coupling. 
Now an increasing number of chemists and physicists are 
taking interest in heterospin ferromagnets^i 2 ^ 



Alternating-spin chains possess elementary excitations 
of dual aspect. In the case of antiferromagnetic coupling, 
the acoustic excitations reduce the ground-state magne- 
tization and are thus of ferromagnetic nature, while the 
optical excitations enhance the ground-state magnetiza- 
tion and are thus of antiferromagnetic nature. In the case 
of ferromagnetic coupling, on the other hand, both ex- 
citations are of ferromagnetic character. Therefore, the 
Schottky-type peak of the specific heat and the minimum 
of the susceptibility-temperature product, which are both 
ferrimagnetic features, are absent from mixed-spin ferro- 
magnets. Nevertheless, mixed-spin ferromagnets and fer- 
rimagnets behave similarly at low temperatures, which 
is the goal of this paper. Employing a quantum Monte 
Carlo method^ 3 , and a modified spin- wave theory^ 2 * we in- 
vestigate thermodynamis of one-dimensional alternating- 
spin Heisenberg ferromagnets with particular emphasis 
on the intrinsic low-temperature properties. 



II. MODIFIED SPIN- WAVE SCHEME 

We consider two kinds of spins S and s (S > s) alter- 
nating on a ring with ferromagnetic exchange coupling 
between nearest neighbors, as described by the Hamilto- 
nian 



N 



Sn-{ 



(2.1) 



Even in one dimension, the conventional spin- wave 
theory^iSSi 2 ^ gives a fine piece of information on the 
ground-state correlationi^&Sii As for the thermal quan- 
tities, however, the low-temperature series expansion 
within the conventional scheme 3 ^ only reproduces the 
leading term of the specific heat and nothing correct for 
the magnetic susceptibility^ Then Takahash: 32 modified 
the spin-wave formalism, imposing a constraint on the 
magnetization, and obtained an excellent description of 
the low-temperature thermodynamics of low-dimensional 



2 



ferromagnets. We develop the modified scheme for 
mixed-spin ferromagnets. 

In order to describe the spin deviation in each sublat- 
tice, bosonic operators are introduced as 



25 - a 1 n a n a n , S* = S - ala n , 



+ = \ 2s - btbn b r , 



S - blb n , 



(2.2) 



where we regard 5 and s as quantities of the same order. 
The bosonic Hamiltonian reads 



H = E 2 +H 1 +H Q + 0(S- 1 ), 



(2.3) 



where E 2 — —2SsJN is the classical ground-state en- 
ergy and Hi is the 0(S l ) quantum correction to it. We 
consider first diagonalizing Tli and then taking TLq into 
calculation perturbationally.24 Via the transformation 

a « = ^E e ~ ifc( "~ 1/4) (4cosfl fc -PlsinO k ), 
b n = -]= J>- ifc (" +1 /4) (at sin 6k + 0l cos 6 k ) , 

(2.4) 

with tan(20fc) = 2 ^/Ss cos(fc/2)/ (5 — s), we obtain 

m=jJ2(^a{a k +u+0tf3 k ). (2.5) 



Here the acoustic (u> k ) and optical {co k ) dispersion rela- 
tions are given by 



uj^ = S + s± y/(S- s) 2 + 45scos 2 (fc/2) 

= S + S ±L0 k , 



(2.6) 



and plotted in Fig. 

Now we proceed to the modified spin-wave scheme in 
an attempt to avoid thermal divergence of the number 
of bosons. At finite temperatures, we replace a\a k and 
PlPk by fi^ = J2n-,n+ n^Pkin" ,n + ), where P k (rT ,n + ) 
is the probability of wT acoustic and n + optical spin 
waves appearing in the fc-momentum state, and minimize 
the up-to-0(5 1 ) free energy 



k <r=± 

+ k sTY, E Pk(n-,n+)\nP k (n-,; 



(2.7) 



with respect to Pk(n , n + )'s under the condition of zero 
magnetization 



0. 



(2.8) 



k ct=± 



together with the trivial constraints 
En-.n+^U"" ™ + ) = 1. Up to CK5 1 ), the mag- 
netic susceptibility and the internal energy at thermal 



equilibrium are expressed as 



3k B T 



k <r=± 



E = E 2 + 

k <y=± 



(2.9) 
(2.10) 



with rit 



[ e (Ju k -»)/k B T _ w here the 5 factors 

of the spins 5 and s are both set equal to g and the 
Lagrange multiplier is determined through the condi- 
tion l|2.8|l . The specific heat is calculated by numerically 
differentiating the internal energy. The perturbational 
correction of O(S ) reads 

(H ) = Tr[H e-^ /feBT ]/Tr[c- Wl / fcBT ] 
JN 



-Pi + Pi 



r 2 



r 2 )r 3 



r 2 vr a 



(2.11) 



with 



r 2 = 



TV ^ 

AT ^ 



5-. 



2V5s 



9 , 

cos -(n fc 



(2.12) 



where we keep /1 unchanged. Indeed fi may be modified 
so as to minimize the up-to-0(5 ) free energy, but such 
a procedure, which is much more complicated, has no 
effect on the low-temperature leading behavior^ 



III. NUMERICAL CALCULATIONS 

Quantum Monte Carlo calculations are presented at 
N = 32 (64 spins). In the case of (5, s) = (1, |), we have 
carried out preliminary calculations at N — 24,32,40. 
Any quantity divided by N does not vary with N be- 
yond its statistical error in the temperature range to 
show. A few million Monte Carlo steps are spent on low- 
temperature calculations, while less than a half million 
steps on high-temperature calculations. The numerical 
precision in the final results is two to three digits. 

In Fig. [3we compare modified-spin-wave and quantum 
Monte Carlo calculations of the magnetic susceptibility. 
They are in excellent agreement over the whole tempera- 
ture range. The analytic calculation reproduces the para- 
magnetic susceptibility [S(S+l)+s(s+l)]N(gn B ) 2 /3k B T 
at high temperatures and reveals the T _2 -diverging be- 
havior at low temperatures, which is later discussed in 
more detail. The susceptibility-temperature product still 
monotonically decreases with increasing temperature in 
contrast with the ferrimagnetic behavior. 34 
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In Fig. [3] we compare the modified spin-wave and 
quantum Monte Carlo calculations of the specific heat. 
The agreement between them is somewhat poorer than 
that found for the susceptibility, but the T 1 / 2 -initial be- 
havior at low temperatures and the spin-dependent peak 
structure at intermediate temperatures are well repro- 
duced by the analytic calculation. The mid-temperature 
structure of the specific heat may be regarded as a func- 
tion of the acoustic excitation band width W~ = 2sJ 
and the optical excitation gap A = 2SJ (see Fig. [Q. 
The heat capacity attributable to the acoustic excita- 
tions and that to the optical excitations may be separa- 
ble when W~ <C A. Observing further calculations in 
Fig. 0] we are convinced that the double-peaked struc- 
ture may appear for S ^ 3s, including practical cases 
(S, s) = (§, i), (2, |), (§, |). Mn(II)Cu(II), Fe(II)Cu(II), 
and Co(II)Cu(II) chain compounds^ have indeed been 
synthesized so far, but they all exhibit antiferromag- 
netic intrachain interaction. The double-peaked struc- 
ture is much more pronounced for ferromagnetic intra- 
chain interactioni^SiS We expect an increased effort to 
design ferromagnetic exchange coupling between alter- 
nating metal ions. 



IV. ANALYTICAL CALCULATIONS 

In order to elucidate the low-temperature thermal be- 
havior, we define the state density function 



(4.1) 



Here we are interested in the gapless acoustic branch and 
expand w~ (x) for small x as 



1 S + s 



(4.2) 



n=0 



A few leading coefficients are given as 

c o = h 

(g - sf + Ss 
Cl ~ ASs(S + s) ' 
_ _ (3S 2 - ASs + 3s 2 ) (5 + s) 2 - 5S 2 s 2 



32S 2 s 2 (S + s) 



(4.3) 



Applying Eq. I|4.2|l and neglecting the optical excitations 



Eq. (EE)) reads 



,1/2 



w^/2Ss(S + s) ^ Q 



/2 



xr(n + - 



n--n)r" 



1 



n — m + — I v 



m=0 



m+1/2 



(4.4) 



where v = —fi/faT, t = k B T/ J , and C( n ) is Riemann's 
zeta function. Solving this equation iteratively, we obtain 
the low-temperature expansion of the Lagrange multi- 
plier as 



,1/2 = 



Co -r(i/2) r t\ 



iry/2Sa(S + a)' 



/2 



c -r(i/2) 



c r(i/2) 



n v / 2Ss{S + s) 

3 



TTy/2Ss{S + T) 



t 3 ^ + 0(t 2 ). 



(4.5) 



Since the magnetic susceptibility and the internal energy 
read 



X^bT 1 S + s 
N(g^i B ) 2 ~ 3^V 2Ss 



/2 



NJ 



xr( n + i)[r(2-„)„-^ 
£^(»-»-5>" 

" oo 

E c ^" +3/2 



m=0 

E-E 2 _ I IS - s 

7T 



(4.6) 



2Ss 



n=0 



xr[n+- 

(-1) 



r(- l --n)v n+1 i 2 



^ ml ^( 



m=0 



n — m H — \v 
2 * 



(4.7) 



the susceptibility and the specific heat are expanded as 



xJ 



1 \i- 2 C(l/2) .--3/ 2 



iV( 3( u B ) 2 

"C(l/2) 



+ 



C 



/2tt 



(S + s) 




(4.8) 



|3C(3/2) fl/2 _ f 



| 4V2^ 



15[(S 2 -g s + s 2 )C(5/2)-4C(l/2)] p/2 



32V2vr 



-OlP), 



(4.9) 



where t = t/Ss(S + s) = k B T/JSs{S + s). 
The O(S ) interactions affect the fourth and higher 
terms and therefore, whether through the Holstcin- 
Primakoff transformation 2 ^ or through the Dyson- 
Maleev transformationjSL 3 ^ we reach the same results 
and (JS2J. 

Numerically solving the thermodynamic Bethe-ansatz 
integral equations for the spin-i ferromagnetic Heisen- 
berg chain, Takahashi and Yamada-i obtained 



xJ 

£(.9Mb) 2 



0.04167t~ 2 + 0.145i" 3/2 + 0.17t~ 1 
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+0(t~ 1/2 ), (4.10) 
= 0.7815t 1/2 - 2.00< + 3.5i 3/2 + Oft 2 ), (4.11) 

Lfc B 

where L is the number of spins. When we set S and s 
both equal to \ , the expressions 1)4. 8|l and (f 4 . 1 C)|> coincide 
in their leading three terms, while Eqs. (|4.9|l and (|4.11|) 
in their leading two terms. The modified spin-wave calcu- 
lations are thus reliable and give rigorous information on 
the low-temperature properties. In the case of arbitrary 
S and s, the leading three terms of Eq. 1)4. 8|) and the 
leading two terms of Eq. (|4.9() coincide with those of the 
STpin-[Ss(S + s)/2] 1 / 3 uniform ferromagnetic chain except 
for a common factor. Considering practical combinations 
of S and s, we may estimate that [S's(S' -I- s)/2} 1 / 3 = 
[1 - (S - s) 2 /{S + s) 2 }^ 3 {S + s)/2 ~(S + s)/2. Thus, 
ferromagnetically coupled alternating spins S and s look 
like a ferromagnetic assembly of virtual spins (S + s)/2 
at low temperatures. 

It is also interesting to compare the expressions ()4.8JI 
and (|4.9|) with those of ferrimagnetic chainsiSi It turns 
out that the spin-(5, s) ferrimagnetic low-temperature 
expansions are obtained by replacing J and s by — J 
and — s, respectively, in Eqs. 14.81) and (|4.9|) . In other 
words, antiferromagnetically coupled alternating spins S 
and s look like a ferromagnetic assembly of virtual spins 
[Ss(S-s)/2} 1 /' 3 = [(5+s) 2 /(5-s) 2 -l] 1 /3(5„ s )/2atlow 
temperatures. The quantity [5's(S'— s)/2] 1 / 3 is less intu- 
itive than the corresponding [S^S 1 + s)/2] 1 / 3 in the case 
of ferromagnetic coupling, but it results in a realistic spin 
quantum number S — s when S is equal to 2s. Equation 
(|4.8|) is nothing but the susceptibility of the spin- (S,s) 
ferrimagnetic chain if we set t for k^T / JSs{S — s) in- 
stead of k-gT I J ~Ss(S ' + s). All in all, the low-temperature 
physics is scaled by S+s in ferromagnetic chains, whereas 
by S — s in ferrimagnetic chains. 



gardless of (S, s) and is qualitatively the same as those 
of uniform ferromagnetic chains. The specific heat qual- 
itatively varies with (S, s) and exhibits a rich structure 
at intermediate temperatures. It may be double-peaked 
for 5* > 3s in general. 

The low-temperature behavior has been revealed an- 
alytically. The thermal quantities are still expanded in 
powers of T 1 / 2 and exhibit ferromagnetic features. The 
conventional spin-wave theory misunderstands the low- 
temperature behavior as series of T. The missing terms 
are reproduced through the modified procedure. Ferro- 
magnetic and ferrimagnetic mixed-spin chains are quali- 
tatively alike at low temperatures. The spin-(S l , s) ferro- 
magnetic chain looks like a ferromagnetic assembly of vir- 
tual spins [Ss{S+s)/2] 1 / 3 = [l-(S-sf/{S+sf] 1 / 3 {S+ 
s)/2 ~ (S+s)/2, while the spin-(S l , s) ferrimagnetic chain 
behaves like that of virtual spins [Ss(S — s)/2] 1 / 3 = 
[{S + s) 2 /(S - s) 2 - l] 1 /3(5- _ sj/2. The present find- 
ings are really complementary to the sophisticated high- 
temperature series-expansion calculations^*^ 

The existent bimetallic chain ferromagnet 
MnNi(N0 9 ) 4 (en) 2 possesses a rather weak intra- 
chain exchange coupling (J//cb — 2K), in which the 
low-temperature thermodynamics revealed here is hard 
to verify. Nevertheless such a pioneering material must 
highly motivate further explorations in both chemical 
and physical fields, as was the case with uniform ferro- 
magnetic chainsi 30 ' 31 i 32 i 38 i 39 ' 40 Besides bimetallic chain 
compounds, several authors 41 made a novel attempt 
to design low-dimensional heterospin systems utilizing 
organic triradicals. Mixed-spin chains contain further 
interesting topics such as dynamic structure factors of 
dual aspect 4 ^ and nuclear spin relaxation through the 
exchange-scattering-enhanced three-magnon process^ 
We hope our study will stimulate further experimental 
investigations into mixed-spin chain compounds. 



V. CONCLUDING REMARKS 

Thermodynamic properties of alternating-spin Heisen- 
berg ferromagnetic chains have been investigated in com- 
parison with heterospin ferrimagnetic and homospin fer- 
romagnetic chains. The magnetic susceptibility is a 
monotonically decreasing function of temperature re- 
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FIG. 1: Single-magnon excitation spectra as the rigorous dispersion relations of the elementary excitations for the spin-(S', s) 
ferromagnetic Heisenberg chains. 
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FIG. 2: Modified-spin-wave (MSW) and quantum Monte Carlo (QMC) calculations of the magnetic susceptibility \ as a 
function of temperature for the spin-(S', s) ferromagnetic Heisenberg chains. 
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FIG. 3: Modified-spin-wave (MSW) and quantum Monte Carlo (QMC) calculations of the specific heat C as a function of 
temperature for the spin-(5', s) ferromagnetic Heisenberg chains. 




FIG. 4: Modified-spin-wave (MSW) calculations of the spe- 
cific heat C as a function of temperature for the spin-(6', s) 
ferromagnetic Heisenberg chains in the cases of S < 3s (a) 
and S > 3s (b). 



